Modulated structures in electroconvection in nematic liquid crystals 
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Motivated by experiments in electroconvection in nematic liquid crystals with homeotropic align- 
ment we study the coupled amplitude equations describing the formation of a stationary roll pattern 
in the presence of a weakly-damped mode that breaks isotropy. The equations can be generalized 
to describe the planarly aligned case if the orienting effect of the boundaries is small, which can 
be achieved by a destabilizing magnetic field. The slow mode represents the in-plane director at 
the center of the cell. The simplest uniform states are normal rolls which may undergo a pitch- 
fork bifurcation to abnormal rolls with a misaligned in-plane director. We present a new class 
of defect-free solutions with spatial modulations perpendicular to the rolls. In a parameter range 
where the zig-zag instability is not relevant these solutions are stable attractors, as observed in 
experiments. We also present two-dimensionally modulated states with and without defects which 
result from the destabilization of the one-dimensionally modulated structures. Finally, for no (or 
very small) damping, and away from the rotationally symmetric case, we find static chevrons made 
up of a periodic arrangement of defect chains (or bands of defects) separating homogeneous regions 
of oblique rolls with very small amplitude. These states may provide a model for a class of poorly 
understood stationary structures observed in various highly-conducting materials ( "prechevrons" or 
"broad domains"). 

PACS numbers: 61.30.Gd, 47.20.Ky, 47.65.+a 



I. INTRODUCTION 

Nematic liquid crystals, the simplest type of intrinsi- 
cally anisotropic fluids, continue to provide model sys- 
tems for a wide variety of interesting nonlinear dynami- 
cal phenomena like optical instabilities Q), flow-induced 
nonlinear waves [0 , critical properties of nonequilibrium 
transitions || , and in particular electrically or thermally 
driven convection instabilities M, Q| (see also ]Bj and ref- 
erences therein). 

In nematics the mean orientation of the rod-like 
molecules is described by the director h. Electroconvec- 
tion (EC) driven by an ac voltage U at frequency ui is 
commonly observed in thin nematic layers sandwiched 
between glass plates with transparent electrodes using 
nematics with positive conductivity anisotropy (p a > 0) 
and negative or slightly positive dielectric anisotropy e a . 
In the well studied planarly aligned case, where n is an- 
chored parallel to the bounding plates along a direction 
which we will take as x (we choose the layer in the x-y 
plane) EC sets in directly from the homogeneous state at 
a critical voltage U c (u>) and leads slightly above thresh- 
old to ordered roll patterns associated with a periodic 
director distortion with the critical wave vector q c {w). 
Here we will only consider the most common situation 
where the bifurcation is supercritical and leads to sta- 
tionary rolls with q c parallel to x (normal rolls (NRs)). 
In the usual low-frequency conduction regime, where the 
wavelength is controlled by the cell thickness, this may 
exclude in particular very low frequencies where the rolls 
at threshold may be oriented obliquely (depending on the 
material). In NRs near threshold the director remains in 
the x-z plane, i.e., perpendicular to the roll axis. 



The investigation of homeotropically oriented cells 
using nematics with manifestly negative dielectric 
anisotropy e a < was initiated rather recently, see 
H @i H eI f° r experimental and jl0| |ll , 12 for theoretical 
work. In this case the director is initially oriented perpen- 
dicular to the layer, i.e., in the z direction, so the system 
is isotropic in the x-y plane. Then the first instability is 
the spatially homogeneous Freedericksz transition where 
the director bends away from the z direction, singling 
out spontaneously a direction c in the x-y plane. After 
the transition the slow, undamped variation of the in- 
plane director c (the Goldstone mode) may be described 
by an angle <j). At higher voltages there is a further tran- 
sition to EC which is in many respects similar to that 
in cells with planarly aligned nematics. However, now 
the Goldstone mode has to be included in the descrip- 
tion even right at threshold. It turns out that the torque 
arising when the in-plane director and the wavevector 
are (slightly) misaligned is destabilizing, i.e., it acts to 
increase the misalignment ("abnormal torque"). Then 
the in-plane director in NRs is not perpendicular to the 
(local) roll axis. These NRs with a misaligned in-plane 
director were termed "abnormal rolls" (ARs) ||. Fur- 
thermore one may expect spatio-temporal disorder right 
at threshold and this has indeed been observed, at least 
in the oblique roll regime, where one expects faster dy- 
namics [p, 11, lj|. For NRs the experimental situation is 
not totally clear JD|. 

By applying an in-plane magnetic field, which now de- 
fines the x direction and exerts an aligning torque on c, 
the disorder at threshold can be suppressed. Indeed the 
situation then is similar to that in the planar case. How- 
ever, for a small field a small abnormal torque will suffice 
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to overcome the magnetic alignment, and then one has a 
transition to ordered ARs, where the in-plane director is 
homogeneously rotated out of the x direction. For NRs 
this symmetry breaking is spontaneous, either to the left 
or to the right, and the transition is described by a su- 
percritical pitchfork bifurcation, which has been verified 
experimentally |TJ|, [TJ] . 

In the planarly aligned case one also has a transition 
to ordered ARs, although this occurs at a distance from 
threshold such that a quantitative description is more dif- 
ficult. Since now the director is aligned at the cell bound- 
aries the distortion of the in-plane director is confined to 
the center part of the cell leading to a twist distortion. 
Incidentally, this is the reason why the phenomenon has 
been identified only recently in planar convection p6[ . 
By applying a magnetic field in the y direction one can 
now destabilize c and thus move the AR transition down- 
ward towards the primary instability. The AR transi- 
tion merges with the primary bifurcation when the field 
reaches the strength of the twist Freedericksz field. When 
the two transitions are near each other a simple reduced 
description can be used. This shows that the two classes 
of systems are in many ways similar. 

Above the AR transition one often observes more com- 
plicated structures with or without defects. In particular, 
in homeotropic EC, modulations of the AR mode have 
been observed which leave the roll pattern (i.e., its phase) 
virtually unchanged |l4|, |l5) . Ideally, such structures can 
be considered quasi one-dimensional (ID) with spatial 
variations only perpendicular to the (normally oriented) 
rolls. We wish to address in particular such structures by 
studying the simplest set of coupled amplitude equations 
capable of describing the AR scenario. These equations 
were first derived for homeotropic systems near thresh- 
old jfl], [TtJ but with a slight generalization they can also 
illustrate the planar case. We here present for the first 
time ID solutions of the appropriate type. We will then 
briefly discuss the destabilization of these ID structures. 
Finally we present a new class of fully ordered 2D so- 
lutions occurring at higher voltage (or smaller magnetic 
fields) involving periodic arrangements of defect chains 
(or bands of defects) . We call them static chevrons since 
they are reminiscent of the dynamic chevrons observed 
in the dielectric range of EC, and more recently, also in 
homeotropic convection [jl3|| . 

Our results could also be of relevance for the higher- 
frequency dielectric regime, where the rolls are very nar- 
row and therefore the orienting effect of planar boundary 
conditions is substantially weaker than in the conduction 
range. Recently the weakly nonlinear description of the 
dielectric regime in planarly aligned systems has been 
investigated in detail |is| . 

In Sec. [n] we discuss the basic equations, which take 
the form of an activator-inhibitor system, and their ho- 
mogeneous solutions. In Sec. [II we study ID modulated 
solutions. Their st ability within the full 2 D equations is 
investigated in Sec. IV A and in Sec. IV B some 2D struc- 
tures, in particular defect-free solutions, are presented. 



In Sec. [v] we present the static chevron states. We con- 
clude by putting our results into a more general context, 
relate them to experiments, and present an outlook. 



II. BASIC EQUATIONS AND HOMOGENEOUS 
SOLUTIONS 



In situations where the in-plane director, described 
by an angle tp measured from the x direction, becomes 
an active mode already near the threshold to NRs one 
can describe the system by the following set of coupled 
Ginzburg-Landau equations for the complex patterning 
mode A and the slowly varying angle tp |ll], [l7t 



fd f A 



+ Q x di+e v (di - 2iq c C l pd y -C 2 ^) 



-g\A\ 2 + wd y tp 



A. 



71 = Kitfip + K 3 8%tp 
+ -[-iq c A*(d e 



Tip 

iq c ip)A + c.c. 



(1) 



(2) 



Here we have chosen the x direction along the wavevec- 
tor of the NRs. The angle ip of the in-plane director is 
measured from the x axis. The validity of Eqs. (|l|,||) is 
restricted to small values of the reduced control parame- 
ter (more precisely e/g <C 1 and angles \ip\ <C 1). Special 
attention should be paid to the sign of the parameter 
T. If T were positive the field (p would be stabilized by 
the roll pattern. However, T is negative, at least for the 
standard nematics which have been used in relevant ex- 
periments Jll], |l7, 18 . This gives rise to the transition to 
abnormal rolls and to interesting dynamical phenomena. 
Note that £, yy tends to zero at the transition from normal 
to oblique rolls at threshold. 

The equations can be justified most convincingly for 
homeotropic orientation near the EC threshold (reduced 
control parameter e <C 1). Overall rotation invariance 
then requires that the three terms multiplying t; 2 y in 
Eq. (Q) combine to (d y — iq c (p) 2 , so that C\ = C% = 1. 
Without the isotropy-breaking term —Tip, which is re- 
alized easily by an in-plane magnetic field (then T = 
X a H 2 ), the angle <p may not saturate. Then one has to 
resort to a globally rotational invariant generalization of 
Eqs. (Q,||) 0,0- We will see that this is not always 
necessary. 

The parameters Ci, C2 were introduced to allow for 
more general situations like in planar alignment (for 
structural stability one needs C2 > Cf > 0). Then the 
magnetic field term in Eq. (^|) actually models the ori- 
enting effect of the boundaries. In the common conduc- 
tion range, where the wavelength is of the order of the 
cell thickness, the orienting effects are sufficiently strong 
so that interesting dynamics of tp sets in at values of e 
which are too large to allow quantitative description by 
Eqs. (0J^)- Then, in particular, singular mean flow has 
to be included g0|. To reduce the damping of p 
one can either apply a destabilizing in-plane magnetic 
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field and/or go to the dielectric range where the effect of 
boundaries is weaker When T becomes too small 

higher-order terms are needed in Eq. (|J) . 

For calculations it is useful to introduce a scaled ver- 
sion of Eqs. (pl||) for e > 



rd t A = [l + d* + d 2 - 2ic 1(l 
d t 4> = B x dl4> + D 2 dU -h<p + 



~cb 2 - \A\ 2 - ivd y <j>}A(3) 



A*(c 2 d y -i<p)A + < 



(4) 



exist in the range < h < 1. They bifurcate from the 
NRs at h= 1 and they break the symmetry <f> — * — cb of 
the equations. For h > the two AR states are global 
attractors with regions of attraction Ao > 0, (j>a > / < 0, 
respectively. For h > 4r/(l + 4r) the AR states represent 
saddles (two real eigenvalues), otherwise spiral points. 

The case h = needs special attention. Clearly the 
whole band 



where 

A = 
x 

r = 
£>i = 

c\ = 



(e/g^A, ^ = e 1 ' 2 ( b/^ yy q c ^C' 2 ), 
t xx £- ll2 x, y = Zyy£- 1/2 y : i = ^ ig t/{e\T\q 2 c ), 
f\T\q 2 c /{2^g), v = O/i&qcy/Qi), 
2K 3 g/(\T\q 2 C x ), D 2 = 2K ig /(\r\q 2 e yy ), 
Ct/y/fy, c 2 = ^C~ 2 , h = 2Tg/(\T\q 2 c s). (5) 



The damping parameter h gives the ratio of the aligning 
(=isotropy-breaking) torque over the abnormal torque of 
NRs. For large values of h one can set cf) = and disre- 
gard Eq. (^). h can be decreased by either decreasing a 
stabilizing magnetic field (increasing a destabilizing field) 
or by increasing e. Below we will show that, keeping the 
aligning torque fixed, one can write 



£ar/£ 



(G) 



where £ar is the reduced control parameter where the 
transition from NRs to ARs takes place. The parameter 
v describes the action of the gradient of the in-plane di- 
rector on the phase of the rolls. Experimentally it can be 
controlled by varying the frequency. In this paper we will 
be concerned with the range v > where NRs are first 
destabilized by the transition to ARs (see below). For 
v < the zig-zag instability comes in earlier. Various 
features of Eqs. ([|[i]) have been analyzed in 11, [l?], |2^| 
and comparison with experiment has given evidence for 
their validity. 

We first discuss the homogenous solutions of (|3|J 
for rolls with modulation wavevector (Q, P) where A — 
A oe i(Qx+Py) ^ Aq > 0, and a constant in space angle 
4> = 4>$. One is left with the dynamical system 



Td t A = [\-Q 2 -P 2 + 2 Cl P4> -4>l-A 2 ]A , 
dtfo = {Al-h)<t> Q -c 2 PAl. 



(7) 



These equations can be classified as an activator-inhibitor 
system with activator A$ (positive linear growth rate for 
not too large wavevector) and inhibitor c/jq. For simplicity 
we will in the following consider Q = (the results are 
easily generalized). For P = (rolls exactly in the x 
direction; we will deal with this case, except in Sec. 0) 
there is the symmetry <fi — > —<p. The basic uniform state 
Ao = 0, 4>o = 0, is an unstable solution of ([?]), and the NR 
state Aq = 1, 0o = O, is stable only for h>l. Finally, the 
abnormal roll (AR) states 



bo = constant, 



(9) 



is solution of the equation. The segment \cbo\ < 1 is re- 
pulsive whereas the regions \4>q\ > 1 are attractive. The 
AR state (for h = 0) separates the two cases. There is 
a separatrix Aq + 1/(1 + t)</)q = 1, which separates tra- 
jectories flowing out of the repulsive segment from those 
coming from infinity. 

The degeneracy for h — is presumably realistic for 
homeotropic EC. It is a consequence of rotational in- 
variance. For planar EC the degeneracy is removed by 
higher-order terms, in particular a term proportional to 
4> 3 in Eq. (Q). Nevertheless it is instructive to study this 
limit where results simplify. 



III. ID MODULATED STRUCTURES 

Next we consider modulated structures that leave the 
roll pattern untouched, i.e., which do not involve the 
phase of the complex field A. This can occur generically 
only for spatial variations along x. Then the equations 
take the form 

rd t A = {dl + 1- A 2 ~cb 2 )A, (10) 
d t cb = [D l d 2 x +A 2 ~h)cb. (11) 

In this section we will study these equations. We will 
choose A > 0. 



A. Single domain walls 

We start by discussing domain walls that connect the 
two variants of AR solutions and their interaction. Near 
the AR transition (h near 1) domain walls attract each 
other, so that modulated states are unstable. This can 
be seen from the fact that for 1 — h <S 1 the amplitude 
can be eliminated adiabatically from Eq. ( |l0| ) leading to 



d t cf>= (DiSP x + l-h-. 



(12) 



.4 



AR 



(8) 



In this equation all modulated states (they can be ex- 
pressed in terms of an elliptic integral) are unstable, al- 
though their lifetimes are exceedingly long due to the 
exponentially weak (attractive) interaction between well 
separated domain walls (see below). 

Surprisingly, for lower damping, the interaction ac- 
quires repulsive parts, so that stable modulated struc- 
tures can emerge. In fact, for h — > the interaction be- 
comes purely repulsive so that only periodic modulations 
exist (see below). 
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FIG. 1: Domain wall solution connecting the two abnormal 
roll states A = h, <j) = ±y/T^H for D y = 0.2 and h = 0.01. 



Actually, for h — a static domain wall solution can be 
found analytically. It reads 

A w {x) = A sech(/3a;), <f>w(x) = </>o tanh(/3x) , (13) 



where 



/? 2 



1 



A 2 



2{l + D 1 )(3 2 



1 + 2D X ' - ^ < 

(14) 

It selects the states A — Aq, (j> — zL(f>o from the continuum 
of states (||). We have checked numerically and could find 
no stable static domain wall solution other than (13). We 
note that since 6n>l (13f) connects stable states. 



We mention that the stationary wall ( |l3| ) is embed- 
ded in a continuous family of moving walls connecting 
inequivalent states (^) . Their study is beyond the scopes 
of the present paper. 

Since for h ^ the walls can connect only to the AR 
states and </>ar — * ±1 {¥" 4>o) f° r h — > 0, this limit has to 
be clarified. An example of a numerically stable (within 
the ID equations) wall for small h is shown in Fig. [y. The 
fields go to their AR values at spatial infinity. However, 
in the wall region the field <j) exhibits an overshoot ap- 



proaching the values ±<f>o of (13). The overshoot becomes 
longer and approaches nearer to ±<fo as h— *0. 

In order to understand this behavior better we study 
the spatial decay of the solution into the AR state. This 
is obtained by linearizing Eqs. ( |To| , [Tl| ) around an AR so- 
lution and calculating the spatial exponents. They are 



P 2 ± = h±Jh 2 



4h(l - h) 



(15) 



For h — > the exponents are complex and tend to zero, 
which explains the slow decay. The exponents p± remain 
complex for damping constants h < h osc = (1 + Z?i/4) . 
We will show that this is the bound up to which domain 
walls may repel and thus form stable bound states. 



B. Interaction of domain walls, modulated 
structures 

In simulations of Eqs. (|l0| , |l"l"| ) one easily finds stable 
stationary solutions with more than one domain wall. In 
particular there are periodic solutions which come in two 
varieties: the ones that preserve the global symmetry 
(f>—> — <fi ("symmetric") and others that do not preserve 
it ( "nonsymmetric" ) . This section is mainly devoted to 
periodic solutions. There are also pulse-type states local- 
ized around one of the AR solutions as well as nonperi- 
odic extended solutions. 

Let us first approach these modulated solutions from 
the side of large separation between domain walls by 
studying their interaction. Repulsive interaction is a pre- 
requisite to form stable periodic states and interaction 
of both signs are necessary to form nonperiodic arrays 
p2| . The problem is formulated as follows: Two walls 
are placed symmetrically around the origin at positions 
±xq. Then the fields have zero derivative at x — 0. We 
now focus on the region x > 0, write the amplitude 
and angle as A — A w (x — Xo(t)) + a(x — xo(t), t), 4> — 
4> w (x — Xo(t)) + ip(x — xp( t), t) and substitute these ex- 
pressions into Eqs. ( JTo| , |lT| ) . In the spirit of small vari- 
ations and slow ( "slaved" ) dynamics we keep the terms 
linear in a and tp, but neglect the ones that are small and 
contain time derivatives 



± tA' w = 

X0<t>' w = 



-{d 2 x 



< w a- {Di dl 



)a 

A n 



2A w (f> w tp 7 
-h)<p, (16) 



where the dot denotes time derivatives and the argument 
of all functions is x — xq. 

To satisfy the boundary conditions one needs a = A w — 



A, 



= —A' w and ip = 



at x = 



0. At the other end x — > oo the perturbations a and <p 
should vanish. Multiplying (]!(]) by the adjoint of the 
translational mode {A' w , —4>' w ) and integrating from zero 
to infinity, we obtain an equation which has on the right 
hand side only the boundary terms at x — 0. On the 
left hand side the integrals may be extended to — oo with 
negligible error leading to 



(< > -t < A' 2 >)x - A' 2 
-Di{4>'w 



+ A W (A W — Aqq) 

+ 4>'M w -4>oo)) (17) 



where < ... >= f°° dx... . This formula can be used 

J — oo 

only when the term in brackets on the left hand side is 
positive, i.e., for 



T < T accel =<<$>/< A® > 



(18) 



Otherwise translation becomes an active mode, i.e., one 
expects spontaneous acceleration of the domain wall. 
This can be shown by retai ning the small time deriva- 
tive terms omitted in Eqs. ([lq). We will here assume 
( |l8| ) to hold (experimentally this appears to be the case) 
and comment briefly on the acceleration instability at the 
end of this section. 
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FIG. 2: A symmetric periodic solution for Di = 0.2 and ft = 
0.1. The dashed line gives the field A and the solid line the 
field <j>. It is found by numerical simulation of Eqs. ( ^o|Jll| ) 
with periodic boundary conditions. 



In the special case ft = and for the wall (|l3|), expres- 
sion ([It]) can be evaluated leading to 



[2(1 + D x ) - tD 1 ]± = 12 D x p e 



-2/3x 



(19) 



Thus, at least for r < r acco i one has Xq > 0, i.e., the 
interaction among walls is repulsive everywhere. This is 
consistent with the numerical observation of stable peri- 
odic solutions and no nonperiodic states. 

For hy^O we may use the asymptotic behavior charac- 
terized by the exponents in Eq. ( |l5| ) to evaluate the right 
hand side of (fl7|). In the range of complex exponents p±, 
i.e., for ft<ft osc , the interaction is repulsive or attractive 
depending on the distance between walls. Thus one ex- 
pects stable periodic solutions in appropriate wavelength 
intervals as well as nonperiodic structures, in agreement 
with simulations. An example of a periodic solution is 
given in Fig. || for a field value h = 0.1. 

Finally, in the monotonic range h osc < ft < 1, the right 
hand side of ( p7| ) is 

n 4 

c2[ Ah{l-h) ~ 1] CX P( 2 P-^")- ( 2 °) 

where c is a factor that cannot be determined from the 
asymptotic analysis and p_ is given in (|l^) . One easily 
sees that the term in square brackets is negative, so that 
the interaction between domain walls is now attractive. 
In this range we expect neither stable periodic nor non- 
periodic solutions, although, as pointed out before, the 
lifetime of modulations may be very long. 

Let us now continue the symmetric periodic solutions 
to small wavelengths. Since they conserve the global 
<p —t —<f) symmetry one expects them to bifurcate from 
the NRs. This is indeed the case. From a simple linear 
stability analysis one finds that NRs are unstable with re- 
spect to periodic modes with wavenumber \q\ < q c where 
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FIG. 3: Stability limits of the symmetric periodic solutions 
in the period (=L) vs ft plane for D\ = 0.2. Below the line 
Li they are unstable against breaking the cj) — > — cj> symmetry, 
leading to ARs. Above the line L2 they are unstable against a 
symmetry-broken periodic solution. The dotted line gives the 
minimum period L m i n = 27r^/_Di/(l — ft) where the periodic 
solutions bifurcate from NRs. On crossing the dashed line 
there is a zig-zag instability which is either short wave (thin 
dashed) or long wave (thick dashed). The dash-dotted line is 
the limit of stability for the ARs. 
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FIG. 4: A non-symmetric periodic solution for Di = 0.2 and 
ft = 0.1. The dashed line gives the field A and the solid line 
the field d>. 



Indeed, the bifurcation is of supercritical type. In Fig. ||, 
which summarizes our results on symmetric periodic so- 
lutions, the minimum period L m j n = 2Tt/q c (h) is shown 
(dotted curve). 

When the (symmetric) periodic solutions are born they 
inherit the instability of the NRs with respect to symme- 
try breaking. For a fixed h and by increasing the period 
we find that they are stabilized at some period L\[h). 
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FIG. 5: A nonperiodic solution for Di = 0.2 and fo = 0.1. The 
dashed line gives the field A and the solid line the field <f>. 



The curve L\{h) can be calculated by a linear stability 
analysis. It diverges to infinity for h — > h osc since above 
this value no stable periodic solutions are expected to ex- 
ist according to the asymptotic analysis given earlier. At 
L\{h) a branch of unstable periodic solutions with bro- 
ken <p — > — <j) symmetry bifurcates from the symmetric 
periodic solutions. We find that the symmetric periodic 
solutions lose stability (for the first time) at L^iji) in 
Fig. | and fr om here on one has a stable nonsymmctric 
periodic solution, where long and short domains alter- 
nate. An example of such a solution is shown in Fig. |] 
for h = 0.1. Clearly this is the result of the nonmono- 
tonic behavior of the interaction between domain walls. 
Presumably the curve L2(h) also diverges at h osc . 

An example of a more complicated nonperiodic solu- 
tion is shown in Fig. |[ Such states have been also ob- 
served in Ref. |22| . There is no definite distance between 
the walls and there seems to be no repeated structure. 
The state can be spatially chaotic. 

Equation Jl7j ) indicates that xq diverges at T = T ajcce \. 
For h = one finds from Eq. ( |l~9| ) 



^acccl — j-. ■ (22) 

Computer simulations show that the static wall is unsta- 
ble, for r > Taccelj to a steadily traveling wall. One sees 
that on crossing the value T aC cei the wall is accelerated 
while the direction of motion is spontaneously chosen 
p3[ . No further study of moving walls will be given in this 
paper. We find numerically that r accc i, as calculated from 
Eq. @, increases with h (r acccl (/i = 0) = 12,r acce i(/i = 
0.1)~16.6,r acccl (/i = 0.5)~26.6 for D X =Q.2). One could 
see this also by inspecting the form of the walls. 



IV. 2D SOLUTIONS 
A. Stability of ID modulations in the plane 

In view of the experimental observation of periodic 
modulations it is essential to study the stability of the 
solutions discussed in the previous section in the context 
of the full Eqs. (||^) in two dimensions. In order to reduce 
the number of parameters we will in this section consider 
the case c\ = c% = 1 appropriate for homeotropic sys- 
tems. Some results on the case when c\ is smaller than 
1 will be given in the next section. 

For the AR solution (j|) one can go through a standard 
linear stability analysis which gives in the long wave- 
length limit that, for i/>0, ARs become unstable against 
zig-zag perturbations at h = 2/3 Jll], [l7| pif . We should 
also recall that ARs are born from NRs as a stable state 
at h = 1 which means that they exist stably for 

h A R = | < h< 1 . (23) 

The value /iar is denoted by a dash-dotted line in Fig. ||. 
We use in this paper parameter values typical for the 
standard substance MBBA. Specifically we take JlTj 

£>i = 0.2, D 2 = 0.5, r = 0.5 v = 0.6 . (24) 

The positive value for v indicates that we deal with fre- 
quencies above the codimension-2 point where the zig-zag 
instability of NRs is not operative |nj [l6[ || || . 

In order to check the stability of the symmetric peri- 
odic solutions found in the previous section against two- 
dimensional perturbations we have performed a numer- 
ical Floquet analysis. Thus we write the perturbations 
a, if of a periodic solution Al(x), 4>l(%) with period L as 

oo 

a(x,y,t) = e at+is v y ^ a n exp{i2nnx / L) , (25) 

n— — oo 
oo 

<p(x,y,t) = e at+ls y y V n exp[i2Trnx/L), (26) 

n— — oo 

where a n ,tp n are constants. We linearize (|]]4|) in the 
perturbations and expand the coefficient functions in the 
linearized equation, which have the periodicity L, also 
in a Fourier series. After truncation a system of linear 
equations is obtained which can be solved numerically 
to give the growth rate cr i (s y ) for the linear mode with 
wave vector s y . The index i runs over the number of 
Fourier modes used. The number of modes necessary to 
achieve a prescribed accuracy is found easily by numerical 
experimentation. 

We have focused on periods which lie between the lim- 
its Li and L 2 in Fig. ^ i.e., to solutions which are stable 
with respect to x perturbations. For the parameter val- 
ues (|24j) our results are included in Fig. |^ (dashed line). 
The thin dashes (for 7.5 13) indicate a short wave 

instability of the corresponding period when crossing the 
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line to the left. The thick dashes (for 5.7<L<7.5) in- 
dicate a long wave instability. It is most interesting that 
the present instabilities occur for a value of h consider- 
ably smaller than /iar which means that there are stable 
periodic solutions for values of h for which the ARs are 
unstable. 

We give the numerical values for the stability limits 
and the wave vector s y of the unstable mode for some 
periodic solutions that we typically use in the simulations 
of this section (they fit into our typical system length of 
64): 

period 8 : 0.23 < h < 0.67 , s v = 0.67 , 
period 9.14 : 0.30 < h < 0.74 , s y = 0.69 , 
period 12.8: 0.42 < h < 0.85 , s y = 0.66. (27) 

The lower limit corresponds to the short wave instability 
along the y direction with wave vector s y . The upper 
limit is due to the instability along the x direction (curve 
LiinFig. §). 

We conclude that by decreasing the parameter h 
(which corresponds to increasing the voltage or decreas- 
ing the magnetic field in a relevant experiment) below 
/iar the ARs become unstable but periodic states with 
certain periods remain stable. They are then expected 
to be observable under appropriate experimental con- 
ditions. As an example we quote the experiments of 
Rcf. [fl4"[ for which the appropriate parameter values are 
close to that we use throughout this section. It is 
reported that for a control parameter corresponding to 
h = 0.4 in our theory, periodic modulations of the di- 
rector orientation are observed. Indeed this value falls 
into the range where ARs are unstable but some periodic 
states are expected to be stable. Concerning the peri- 
odicity, the tendency towards shorter wavelengths with 
increasing voltage reported in |l4"| is consistent with our 
results. 



B. 2D modulated structures 

We now proceed to investigate how the system evolves 
once the periodic states are destabilized at small values 
of h, that is, once we cross to the left of the dashed 
line in Fig. |3[ For this purpose we have used a pseu- 
dospectral algorithm which simulates the time evolution 
of Eqs. (^Q) which we describe briefly in the following. 
The linear part of the equation is tranformed into Fourier 
space and the analytic formula for its time evolution is 
implemented in the algorithm. For the nonlinear part we 
work in real space and integrate it in time using a vari- 
ation of the Eulcr method. We use periodic boundary 
conditions in both space dimensions which are practi- 
cally enforced by the use of Fourier modes. We typically 
use 128x128 modes in the two space directions while our 
physical space has dimensions 64 x 64 units. Due to the 
special treatment of the linear part, pseudospectral algo- 
rithms allow for a relatively large time step. We typically 
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FIG. 6: A two-dimensional periodically modulated pattern 
without defects from a simulation of Eqs. ([],§) for the param- 
eter values (p^), C\=C2 — 1, and /i = 0.28. Initial conditions: 
a slightly perturbed periodic modulation state with period 
9.14 . Shown is a grey-scale representation of the field \A\. 
White corresponds to the maximum value (\A\ = 0.85) and 
black corresponds to the minimum =0.12). The physical 
dimensions of the space are 64x64 and we have used periodic 
boundary conditions. 




FIG. 7: The field (f> for the state of Fig. ^j. White corresponds 
to the maximum value and black corresponds to the minimum 
value (4> m = ±0.998). 



take St — 0.05 when we use the parameter values (|24|). 

We start a simulation of Eqs. (||,(l|) with a state that 
is periodic in x, e.g. with period 9.14. The state is sta- 
ble in the parameter range given in Eq. (|27|). When we 
reduce the field h slightly below the value 0.30 we find 
that the pattern is destabilized and modulations along 
the y direction appear. The pattern becomes periodic in 
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FIG. 8: A reasonably developed defect lattice from simula- 
tions of Eqs. (|,^) for the parameter values ci =C2 = 1, 
and ft = 0.27. Initial conditions: a slightly perturbed periodic 
modulation state with period 9.14. Shown is a grey-scale rep- 
resentation of |j4|. White corresponds to the maximum value 
and black corresponds to the minimum (\A\ = 0). The state 
is dynamic but always remains close to the one shown. The 
physical dimensions of our space are 64x64 and we have used 
periodic boundary conditions. 



both spatial dimensions. We give an example of such a 
state in Figs. ^, ^ where the fields \A\ and </> are shown, 
respectively. The state describes modulations in 2D and 
has no defects. It was obtained for h = 0.28. We note 
that the state in Figs. is not fully static. It persist 
for a long time but it is eventually modified through the 
creation of defects. The system presents persistent dy- 
namics until the end of our simulation, nevertheless, the 
2D correlations are preserved to a large extent. 

In general, the 2D modulated states reached slightly 
below the destabilization of ID periodic modulations are 
delicate and we have not been able to find a truly static 
one. The near-periodicity of the final states in both spa- 
tial directions is readily seen in the figures resulting from 
our simulations. 

An interesting point is that the defects which are typi- 
cally created in these processes also show strong 2D cor- 
relations. In fact a slightly disordered defect lattice is the 
usual outcome of such numerical simulations. An exam- 
ple of a reasonably developed defect lattice is presented 
in Fig. H through the field \A\. It was obtained by start- 
ing the simulation with a periodic state with period 9.14. 
We used the value h = 0.27 and the other parameters as 
in The initial ID periodic state is then unstable 

and evolves into the defect lattice. The time evolution 
shows that the state has some substantial dynamics but 
it always remains close to the picture of Fig. ||. Finally, 
in Fig. g we give the <f> field of the defect lattice of Fig. [|. 
The <fr field is predominantly varying in the x direction 




FIG. 9: The field <j> for the state of Fig. g. White corresponds 
to the maximum value and black corresponds to the minimum 
value. 




i W * 4 , * * * * * * 1 



FIG. 10: A disordered state: grey-scale representation of the 
field \A\. The parameter values used are ci=C2 = l and 
h = 0.1. The space dimension are 64x64 and we have used 
periodic boundary conditions. 



with small modulations in y. 

On further reducing the parameter h the resulting pat- 
tern becomes progressively more disordered. A defect 
chaotic state is eventually obtained for small values of the 
parameter. In Fig. 



10| we show such a state for h = 0.1. 
These defect-chaotic states deserve to be studied on their 
own right but a full investigation of this problem is be- 
yond the scopes of the present paper. 

It is not easy to predict in what way the states that we 
study here would appear in an experiment. Details of the 
experimental procedure could be important and the final 
result can be complicated as is indicated, e.g., by Fig. 7 



9 



in H] where roughly periodic states without defects are 
found. In order to observe defect lattices in homeotropic 
cells one should remain at small values of e so that the 
director angle does not become very large. Otherwise 
a transition to CRAZY rolls occurs first, which entails 
formation of disclinations in the director field JIJ, [L5| . 
Also, at higher values of e, mean- flow effects not included 
in our simple description become important. 

Let us describe the protocol of a numerical simulation 
which may give some guidance for experiments. The sim- 
ulations are performed on a 64x64 domain. We start at 
a field h above /iar (h = 0.7) with small spatially ran- 
dom initial conditions of the fields A, (f>. Experimentally 
this corresponds to jumping from below the EC thresh- 
old directly into the region of stable ARs. Typically, 
large domains with both AR states appear, separated by 
walls in the x direction. We then jump to h = 0A where 
ARs are destabilized while defects are created and a state 
with 2D correlations sets in. The final state is roughly 
periodic with period 12.8 in the x direction (5 periods 
in the domain) and 16 in the y direction (4 periods in 
the domain). The defects are also roughly ordered in a 
lattice. After this we increase h gradually in steps of 
<5/i ~ 0.02 and let the system relax in every step for 1000 
time units. Although the system never relaxes to a static 
state, we observe clearly correlations in both spatial di- 
rections with the defects approximately ordered in arrays 
until h = 0AA is reached. For h>0A6 the defects annihi- 
late and the system relaxes to a static ID periodic state 
with period 12.8. The process described above seems ro- 
bust in our simulations. Note the interesting fact that 
there is a range of stable coexistence of the fully ordered 
ID state (it is stable down to h = 0.42) and the 2D so- 
lution. 

Defect lattices have in fact been observed in planar 
EC (2^, |27j at fairly high frequencies. A theoretical de- 
scription for these systems by a much more elaborate 
quantitative theory has been developed recently |2(|. 



V. STATIC CHEVRONS 

Here we will relax the condition c\ = C2 = 1 by reduc- 
ing ci (slightly) below 1. This has no influence on the ID 
solutions and their stability with respect to x dependent 
fluctuations and increases the range of stability with re- 
spect to y variations. In fact the critical value of h for the 
zig-zag destabilization of the periodic solutions (dashed 
line in Fig. ^) is reduced. For example, the period 9.14 
is stable in the range 



0.19 <h< 0.74, when ci = 0.9, c 2 = 1 . 



(28) 



which should be compared to Eq. (|27j). For values of h be- 
yond the instability the states that we obtain are similar 
to the ones found above for ci = 1. Namely, we observe 
states with 2D order close to the instability point (with 
some hysteresis), which get progressively less ordered as 
h is decreased. When c\ is decreased further the stable 
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FIG. 11: Static chevron state: grev-scale representation of the 
field \A\. The parameter values used are (|24|), ci =0.9, C2 = 1, 
and h = 0. The physical dimensions is 64 x 64 and we have 
used periodic boundary conditions. 




FIG. 12: The field for the static chevron state of Fig. 0. 



range of the periodic solutions extends to progressively 
lower values of h. In fact, for c\ < 0.77 the ID solutions 
with period 9.14 are stable down to h = 0. 

For h = we obtain a novel static pattern containing 
lines (or bands, i.e., multiple lines) of defects along the y 
direction which resemble the usual chevron states, except 
that those are dynamic. In Figs, [ll] and |l2] such a solution 
is presented for c\ — 0.9. The defects (zeros of \A\) are 
located at the dark points in Fig. O. Within one band 
all defects have the same topological charge and this is 
alternating from band to band. The dark regions between 
defect bands imply a small value of \A\. As one sees from 
Fig. 12 the field varies periodically essentially only in x, 
and thus is very reminiscent of the ID patterns discussed 
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before. We call these new patterns "static chevrons". 
They are formed from random initial conditions. 

The static chevrons persist for small nonzero values of 
the parameter h. If we let the chevron states be created 
for h = and then increase h the static chevrons persist 
up to h< 0.015 for the parameters used here. Increasing 
h further (h = 0.02) we observe that more defect bands 
are created, that is, we have chevrons with a shorter pe- 
riod, but these are now the dynamic chevrons presented 
previously for c\ = 1 0, |28|. For even larger values 
of h one reaches the much less correlated defect chaotic 
state discussed before. For c\ closer to 1 the correlations 
become weaker. 

To gain an understanding of static chevrons we first 
note that the states between the defect bands should be 
interpreted as (approximately) homogeneous states with 
nonzero wavenumber P in the y direction, i.e., with com- 
plex amplitude 



A = A e ±lPy . 



where the sign in the exponent alternates between neigh- 
boring regions and the magnitude of P is related to the 
number of defects p per unit length in one band through 
the relation 



P 

P = — ■ 

7T 



(30) 



This follows directly from the fact that each defect con- 
tributes a phase change of ±2-7r, depending on its topo- 
logical charge. In the simulations we find that and P 
take values close to 



±- 



1 



and Pa = c\i 



(31) 



in the region between the defect bands while Aq is small 
as was already mentioned. 

The uniform solutions are described by equating the 
right-hand sides of Eqs. (Jfy to zero. The resulting cubic 
equation for O can be written as 

(0o - c 2 Pf - 200(1 - Cl c 2 )(0 o - c 2 Pf 

+ [(1 + cj - 2 Cl c 2 )<j>l ~ c 2 2 }(cb - c 2 P) + hc\ O = .(32) 

In the limit h — > one of the solutions tends to 
1 + c| — 2c\c 2 2 



Al = l- 



and the other two tend to 



» , 



90 



Pc 2 



A Q = 0, cf> =c 1 P±Jl-P 2 (l-c 2 1 ) 



(33) 



(34) 



Here one has to impose the additional restriction that 
0o / (0o — c 2 P) > 0, since otherwise Aq is not positive for h 
slightly away from zero. In Fig. |l3| is sketched for these 
limiting solutions as a function of P (solid lines). The 
solution (|33|) is the continuation of the NRs to nonzero 
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FIG. 13: The solutions of Eq. (|32|) (with A% >0) in the (P - 
(29) 0o ) plane for two values of the field h = (solid line) and h = 
0.1 (dashed line). The parameter values are c\ — 0.9, C2 = 1. 
The point A has coordinates 0a ~1/ \J 1 — cf , Pa = ci0a and 
the point B has b =02. /yl + c\ — 2ciC2, Pb = 4>b / 'o*. 



P. Similarly, ( p4|) represents the continuation of ARs. 
The two lines join at point B in Fig. [l^. 

We can now see that the states between the defect 
bands given in Eq. (|3l| ) correspond (approximately) to 
the continued ARs of Eq. ( [34] ) with the maximal \<f>o\. 
Thus the defect bands connect between the two symme- 
try degenerate versions of this state. A simple stability 
analysis explains why this state is selected. We write 
A = Aq + a and = 0o + an d linearize Eqs. (^||) in 
a, ip. For Aq = the two equations separate and we 
obtain the following expression for the growth rate of the 
mode involving a 



TO 



1 



2cic 2 4> s y 



c 2 (, 



(35) 



One easily checks that all states are unstable against 
fluctuations along y except for the state (|3l|), which is 
marginally stable. The latter is denoted by the letter A 
in Fig. H 

The static chevron states which are periodic and de- 
scribe spatial oscillations between the states ( |3l| ) need 
not be themselves marginally stable. On the contrary, 
they are numerically found to be quite robust. This ob- 



servation is in agreement with the results of Sec. IV A 
that the ID periodic states are stable even for parameter 
values for which the uniform states are unstable. 

The effect of the value of ci on the chevron states is 
rather profound. This is seen by the (approximate) rela- 
tion giving the density of defects in the chevrons which 
can be found using Eqs. ( |30| ) and (|3l"|): 



P = 



ci 



ttv/t 



(36) 
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Thus, as ci approaches 1 the defect bands become 
broader (the separation between defects inside a band is 
rather independent of parameters). In the limit c\ —* 1 
one expects existence of an infinite lattice of defects with 
the same polarity, but this is not accessible numerically, 
and presumably also not experimentally. On the other 
hand, as c\ decreases, the number of defects per unit 
length decreases quickly reaching a situation with a single 
chain in a band. Decreasing c\ further one reaches a criti- 
cal value below which no defects are created. Instead one 
has the ID periodic modulations without defects which 
are here stable also for h = (see above). 

We have mentioned already that static chevrons persist 
to slightly nonzero h. In order to understand the defect- 
free regions for this case we have studied the homoge- 
neous solutions as obtained from Eq. (Q) . In Fig. [l3] we 
also show 4>o for the case h = 0.1 (dashed). Analyzing the 
stability we found that there is a small interval around 
the states with maximal \4>q\ which are stable. This shows 
that the disappearance of static chevrons for increasing 
h is not a result of the homogeneous regions becoming 
unstable, but rather the defect bands destabilize, which 
is also observed in simulations. 



VI. CONCLUDING REMARKS 

Motivated in particular by experiments in electrocon- 
vection in homeotropically aligned nematics H, [ll| we 
have studied some classes of modulated solutions of the 
Ginzburg-Landau equation for the complex order param- 
eter A describing the bifurcation to a stationary roll pat- 
tern, coupled to a weakly damped (or even undamped) 
homogeneous mode 4> describing the orientation of the 
in-plane director, see Eqs. (|l],||) (unsealed) or Eqs. (d||) 
(scaled). The most important parameter h in these equa- 
tions gives the ratio of the aligning torque on the director 
over the destabilizing torque of normal rolls. The latter 
is proportional to the supercriticality parameter s. An- 
other important parameter v in the equations character- 
izes the action of the gradient of the in-plane director on 
the phase of the rolls, which can be controlled experi- 
mentally by varying the frequency. Our study is relevant 
for v > 0, which is the case in the upper- frequency con- 
duction range and presumably also in the dielectric range 
& 

The solutions we considered are characterized by mod- 
ulations in the direction perpendicular to normal rolls 
(parallel to the wave vector). There exists a surpris- 
ingly rich spectrum of stable solutions of this type even in 
the range of control parameter where most homogeneous 
states have lost stability. The simplest type has only vari- 
ations in the direction of the wavevector (ID structures). 
When such variations arise they lead to the creation of 
defects, except in a small control parameter region where 
2D defect-free modulations may exist metastably. The 
solutions with defects range from essentially fully spa- 
tially ordered defect lattices to defect chaos with various 



degrees of spatial correlations. Their detailed study ap- 
pears interesting, but is beyond the scope of this paper. 
The complexity of the solutions, as well as their dynam- 
ics, generally increase with decreasing h. For very small 
values of h, however, defects show the tendency to or- 
der along chains (or bands) yielding chevron structures. 
One mechanism for the generation of dynamic chevrons 

~1 We 



has been related to a Turing instability |17 



have found for c\ < 1 a new type of fully ordered, static 
chevrons that appear more related to the ID structures. 

Focusing on the in-plane (= c) director field, the mod- 
ulations occur along the direction of prealignment of the 
in-plane director (from this point of view they give rise 
to "bend" distortions). In fact, there is a long history of 
observations of such modulated states under an electric 
field in highly-doped_MBBA g [§, [H], |3§ and other 
nematic materials 



33| , p4| at high frequency, often with- 
out detectable convection rolls. They come under the 
name of "wide domains" [3(], |33| or "prechevrons" 
fjl| , |32| . Since such structures cannot be explained purely 
statically, it seems reasonable to assume that they are 
secondary structures of the type discussed here. The 
mechanism for the generation of the primary roll pat- 
tern, which may not be visible, is presumably different 
from the usual Carr-Helfrich mechanism. In fact, it was 
shown that the driving mechanism persisted above the 
nematic- isotropic transition |52|. Moreover, in [[33| it 
was shown that a treatment of the bounding plates by 
tensides led to a considerable increase of the frequency 
range where the wide domains appeared. These findings, 
which are consistent with some old results jj5|, indicate 
that an isotropic mechanism with rolls confined to the 
boundaries is involved. 

It was found that even after the wide domains have 
formed the usual electroconvection would set in at the 
expected threshold ||l], |3^] . In the dielectric regime this 
would lead directly to perfectly ordered chevrons that 
have an appearance like the static chevrons presented 
above. 

While states modulated only along the direction of the 
wavevector have been observed to some extent, the ob- 
servation of corresponding two-dimensional modulations 
without defects is lacking. This may be due to the small 
parameter range where such 2D modulations occur and 
also because of their apparent metastability. An exper- 
imental effort in this direction appears interesting. The 
present results suggest that the objective should be to 
observe successively the ID and 2D modulations and the 
chevron states and to test if they can be related to each 
other according to the present theory. 

Equations (§,§) also allow for solutions which are mod- 
ulated only along the rolls (perpendicular to the prealign- 
ment of the in-plane director, giving rise to a "splay" 
distortion). Then the phase of the patterning mode (i.e., 
the phase of A) comes into play, and the periodic solu- 
tions in fact correspond to zig-zag patterns. Such states 
are of relevance for smaller values of v, in particular for 
v < 0. Near the zig-zag instability of normal rolls, which 
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occurs at h = 1 — v, and thus precedes the transition to 
abnormal rolls for v < 0, the periodic patterns are unsta- 
ble. For smaller values of h (larger e) there are various 
types of periodic patterns. The scenario is enriched by 
the fact that the usual defects (point dislocations) can 
extend in the direction perpendicular to the rolls to form 
phase slip lines [|TJ . 

In order to extend the treatment to larger e one has 
to use more complicated equations derivable from an 
extended nonlinear analysis, which includes in particu- 
lar mean flow. For planar convection in the conduction 
regime the analysis has been carried out and relevant 
solutions involving in particular realistic defect lattices 
have been studied pO| . For the dielectric regime an 
extended weakly nonlinear study was presented in [jl8| . 
The system appears interesting in view of its parameters 
ci = C2 = v = 0.98 calculated under neglect of flexoelec- 
tric effects, which are of more relevance in the dielectric 
regime (note that in Ref. Jl3| the Cj correspond to our 
Ci). There is hope to find static chevrons if the align- 



ing torque could be made sufficiently small by applying 
a destabilizing magnetic field. One might use this tech- 
nique also in the conduction range where c\ is substan- 
tially smaller than 1 E^j. However, one may then have 
to include higher-order terms in Eq. fl). 
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